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The WISE view of the disc-torus connection in z ^0.6 Active 
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ABSTRACT 

We selected all radio-quiet AGN in the latest release of the Sloan digital sky survey quasar 
catalog, with redshift in the range 0.56-0.73. About 4000 (~80%) of these have been detected 
in all four IR-bands of WISE (Wide-field Infrared Survey Explorer). This is the largest sam- 
ple suitable to study the disc-torus connection. We find that the torus reprocesses on average 
~ 1/3-1/2 of the accretion disc luminosity. 
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1 INTRODUCTION 



Since the observations of NGC 1068 in polarized light by An- 
tonucci & Miller (1985), showing the presence of broad permitted 
lines in emission, the idea of the unification scheme of radio-quiet 
Seyfert galaxies and quasars emerged (for reviews see Antonucci 
1993; Robson 1996; Peterson 1997; Wills 1999; Krolik 1999). The 
simplest version of the scheme assumes the presence of a dusty 
"torus" surrounding the central regions of the Active Galactic Nu- 
cleus (AGN) intercepting a fraction of the illuminating accretion 
disc radiation and re-emitting it in the infrared. If the absorption 
is due to dust, there is a natural temperature scale in the system, 
since dust sublimates for temperatures greater than ~ 1500 K, cor- 
responding to a peak in the corresponding black body spectrum at 
Up = 3.93 kT/h ~ 1.2 x 10 14 Hz (or A p ~ 2^m; the 3.93 factor 
is appropriate for the peak in the vL v spectrum). The torus origin, 
stability, structure (see e.g. Krolik & Begelman 1988) and its very 
presence in both highly luminous radio-quiet quasars and in low 
luminosity radio loud sources is under debate. The amount of re- 
processed IR radiation seems to become smaller for larger optical 
luminosity in radio-quiet objects (i.e. "receding torus", Lawrence 
1991), while, for radio sources, the absence of broad emission lines 
in low power FR I radio-galaxies (and BL Lacs) could be intrinsic, 
and not due to an obscuring torus (Chiaberge, Capetti & Celotti 
1999). Along the years, the idea of a simple and uniform "dough- 
nut" around the accretion disc has been replaced by a clumped ma- 
terial, possibly outflowing (or inflowing), as envisaged and mod- 
eled by many authors (see e.g. Elvis 2000; Risaliti, Elvis & Nicas- 
tro 2002; Elitzur & Shlosman 2006; Nenkova 2008). 

The existence of the unifying scenario based upon intrinsi- 
cally equal but observationally different AGN is also at the base of 



synthesis models of the X-ray background (Setti & Woltjer 1989; 
Madau, Ghisellini & Fabian 1994; Comastri et al. 1995; Gilli, Co- 
mastri & Hasinger 2007), since also the X-rays are partly absorbed, 
and partly (Compton) reflected by the torus (Ghisellini, Haardt & 
Matt 1994): for large viewing angles, the observed X-ray emission 
becomes very hard, as required to fit the X-ray background. 

The covering factor of the absorbing material forming the 
"torus" is not well known. Estimates come from direct observa- 
tions of optical and IR AGN, as well from statistical considerations 
concerning the number of type 1 and type 2 AGN. In the first case, 
the studies were hampered up to now by the relatively small sam- 
ples of objects (especially in the IR) suitable for a combined study 
(see e.g. Landt et al. 201 1 for a sample of 23 objects observed spec- 
troscopically in the optical and in the IR, down to ~ 3/im). 

In order to study the accretion disc-torus connection in AGN, 
we need to collect the largest group of radio-quiet AGN with re- 
liable detections of the IR luminosity and an optical spectrum to 
characterize the accretion disc features. The Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000) and the Wide-field Infrared Sur- 
vey Explorer (WISE; Wright et al. 2010) are the catalogs with the 
widest number of objects in these two bands, hence they are the 
most appropriate for our study. WISE provided photometric obser- 
vations in 4 IR bands (3.4, 4.6, 12 and 22 /im) for half a billion 
sources (all sky) with fluxes larger than 0.08, 0.11, 1 and 6 mly in 
unconfused regions on the ecliptic in the four bands. The sensitiv- 
ity improves toward the ecliptic poles due to denser coverage and 
lower zodiacal background 

We adopt a flat cosmology with Ho — 71 km s _1 Mpc -1 and 
n M = 0.27. 
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2 SAMPLE SELECTION 

We consider the fifth edition of the SDSS Quasar Catalog (Schnei- 
der et al. 2010), containing 105,783 quasars with magnitude 
smaller than M v band = -22 (i.e. i/L„(5100A) ~ 10 44 ergs -1 ), 
at least an emission line with FWHM> 1000 km s -1 and a reli- 
able spectroscopical redshift. Continuum and line luminosities (as 
well as many other spectral properties) in SDSS spectra have been 
measured by Shen et al. (2011; hereafter Sll). We will use these 
data to estimate optical bolometric luminosity of the sources, fol- 
lowing two independent methods. The first one relies on the 3000A 
continuum luminosity: L^oi = 5.16 x i/L„(3000A) (e.g. Elvis et 
al. 1994; Richards et al. 2006). The second one will be used for a 
consistency check of our results, and relies on H/3 and Mgll line 
luminosities ($37T). In the following we will assume that the bolo- 
metric luminosity equals the accretion disc luminosity. The super- 
script "iso" reminds that they are derived under the assumption of 
isotropic emission. 

The requirement that all sources are observed in the rest frame 
2500-5500 A range (to comprise both the H/3 and Mgll lines and 
the continuum at 3000 A) sets our first selection criterion. Given 
the wavelength coverage of the SDSS, we require a correspond- 
ing redshift range: 0.56 < z < 0.73. The Sll catalog has also 
been cross-correlated with the Faint Images of the Radio Sky at 
Twenty-centimeter survey (FIRST; Becker et al. 1995) and hence 
Sll include in their sample the radio fluxes. The flux limit of the 
FIRST sample is ~1 mjy at 1.4 GHz. Therefore, we can select the 
radio-quiet quasars as those objects observed by the FIRST with- 
out a detectable radio flux. The radio-quiet requirement ensures the 
absence of a contamination from the jet in the wavelength intervals 
of interest. After the radio-quietness and the redshift selections, 
we are left with 5122 sources. We have cross-correlated this sam- 
ple with the WISE All-Sky source catalog requiring that the opti- 
cal and IR positions are closer than 2 arcsec (5082 sources), and 
selecting only those objects with detections in all the four WISE 
IR-bands, to have the most complete IR luminosity information. 
This last selection leaves us with a sample of 3965 WISE-detected, 
radio-quiet type 1 AGN in a redshift range z=0.56-0.73. 



Table 1. Center wavelengths and frequencies of the four WISE 
bands, and value of Am needed to transform the magnitudes 
given in the Vega system to the AB one. See Cutri et al. (2012) in 

http://wise2. ipac .scaltech.edu /docs /release/alls ky 



band 


A [fim] 


log Freq [Hz] 


Am 


1 


3.435 


13.35 


2.699 


2 


4.6 


13.63 


3.339 


3 


11.56 


14.03 


5.174 


4 


22.08 


14.16 


6.620 



Table 2. Mean and standard deviation of IR luminosities and spectral slopes 
in the four WISE bands for the whole sample, and the three subsample 
described in [[3] Luminosities are in units of erg s — 1 , 





Band 


log vL v 


a 


Whole 


1 


44.87±0.26 


1.3±0.5 


sample 


2 


44.92±0.29 


0.9±0.3 




3 


44.87±0.29 


1.4±0.5 




4 


44.99±0.27 




Sub A 


1 


44.88±0.20 


1.2±0.4 




2 


44.74±0.17 


0.9±0.3 




3 


44.77±0.15 


1.5±0.4 




4 


44.74±0.15 




SubB 


1 


45.01±0.20 


1.4±0.4 




2 


44.91±0.19 


0.9±0.2 




3 


44.96±0.16 


1.4±0.4 




4 


44.91±0.15 




SubC 


1 


45.15±0.22 


1.6±0.3 




2 


45.09±0.17 


0.8±0.2 




3 


45.16±0.15 


1.2±0.4 




4 


45.09±0.13 





3 DATA ANALYSIS AND RESULTS 

For all the 3965 sources in our sample we computed the IR flux 
in the four WISE bands by first transforming the observed (Vega) 
magnitudes in the AB systems setting tuab — m + Am, with Am 
given in Tab. Q] In the AB system, the flux-magnitude relation is 
simply: 



logF = - 



m A B + 48.6 
2.5 



(1) 



where the flux density is measured in erg s 1 cm" " 2 Hz -1 . The 
integrated IR luminosity is computed by assuming a power law 
spectrum between two contiguous bands, and summing the con- 
tributions in all the three intervals. The slopes of the power laws 
are given by: 



m-AB.i — mAB.j 



2.51og(A i+1 /Ai) 
The integrated luminosity in each interval is 
v\L v . 



(2) 



(3) 



Finally, the integrated luminosity is L}^ = L2.1 + + L4.3. As 
discussed in $2] the bolometric luminosity is computed as L^oi = 
5.16x^L„(3000A). Again, the "iso" superscript reminds that these 
quantities are computed assuming isotropic emission. 

The ratio R — Lf^/L 1 ^ is approximately constant (~ 0.3, 
Tab. [5] Fig. El) and will be used in Sj4]to estimate the torus covering 
factor. The bolometric and IR luminosities of all sources show a 
well defined correlation over at least 1.5 dex, as shown in Fig.Q] 
We performed two least squares fits by taking at first x — log 
and y — log Lf£, then inverting the variables. We took the bisector 
as the best description of the correlation: log LJr oc 0.83 log L™°x. 
The slope, being smaller than one, suggests that IR luminosities 
become smaller at larger optical luminosity (receding torus). Simi- 
lar results have been found using independent methods by e.g. Ar- 
shakian (2005) and Simpson (2005). 

To provide a deeper insight on the disc-torus connection we 
select three subsamples according to L^°f, we will refer to these 
subsamples with letters A, B, C. Tab. |2]lists the vL v IR luminosi- 
ties in the four WISE bands for the whole sample and for the sub- 
samples A, B, C, together with the average spectral indices. Tab. [3] 
reports the average and the standard deviation of L^oi and Lf£ to- 
gether with their ratio R — LlfJ /L^°i for the whole sample and for 
the A, B, C subsamples. For the latter, instead of the standard devia- 
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Figure 1. Comparison of bolometric luminosity LV°, and integrated IR lu- 
minosity LJ|^ as measured by WISE. Contour levels are at 10%, 50%, 68% 
and 95% of total source number (3965). The dotted line identifies equal lu- 
minosity. The solid line (log LJ^" <x 0.83 log is the bisector of the 
two least squares fitting (dashed) lines. 




Figure 2. Histogram of R, as computed using integrated IR luminosity Lj^? 
and both the continuum-based (red solid line) and BLR-based (green solid 
line) bolometric luminosities. Dashed lines show the cumulative fraction 
(values on right axis). 



tion of I/bob we gi ve the logarithmic width of the considered lumi- 
nosity bin. Note that sources in these subsamples account for only 
~l/3 of the entire sample. Dropping 2/3 of the sample was neces- 
sary to significantly separate the bolometric luminosity classes. 

3.1 Consistency with broad emission lines 

The estimates given above do not take into account that the ob- 
served optical continuum can include different components, be- 
sides the disc emission. As a consistency check, we use an alter- 
native method to derive the disc (bolometric) luminosity, by using 
the luminosities of the H/3 and the Mgll broad lines, always present 
in the SDSS spectra in our redshift selection. For radiatively effi- 
cient discs, indeed, the overall luminosity of the broad line region 
(BLR), I/blr, is a proxy of the disc luminosity Lboi, since on av- 
erage Lblr ~ T^bob where the factor 7 is directly connected to 
the BLR covering factor (see e.g. Baldwin & Netzer 1978; Smith 
et al. 1981). In turn, Ld should be equal to Lboi (real, not isotrop- 
ically equivalent). Estimates of 7 lie in quite large ranges, histor- 
ically between 0.002 and 0.35 (according to Baldwin & Netzer), 
but preferentially < 0.15 (Smith et al. 1981). Typically, an aver- 
age value 7 ~ 0.05 — 0.1 is assumed. Lblr can be calculated 
from individual broad line luminosities, as in Celotti et al. (1997). 
Specifically, setting the Lytv flux contribution to 100, the relative 
weights of the Ha, H/9, Mgll and CIV broad lines are 77, 22, 34 
and 63, respectively (Francis et al. 1991). The total broad line flux 
is fixed at 555.8. Since all our sources have measured H/3 and Mgll 
line luminosity, we average for each object the estimates of Lblr: 



Lblr = — 



555.8 _ , _ _ „, 555.8 _ , , , TT v 
L ^P) + -T7- L(MgII) 



22 



34 



(4) 



A good agreement between the continuum-based and the 
BLR-based bolometric luminosities is obtained using L b s °j = 
Lblr/0.041. This value is found for sources in the highest lumi- 
nosity subsample (C), for which we do not have spurious contribu- 
tions from components other than the disc (e.g. host galaxy). Fig. 
[2] shows the histogram of the ratio R for all sources (green solid 
line). This histogram not only has the same average of the distri- 
bution of R based on the 3000 A luminosity (which is expected, 
given our assumptions), but the two distributions are similar for all 



R, implying that broad lines are a good proxy to compute the total 
disc luminosity, and that our estimates of R are reliable. 



4 THE COVERING FACTOR OF THE TORUS 

Consider the simplest case of a doughnut-shaped torus with open- 
ing angle 9t, as measured from the symmetry axis. The covering 
factor c is defined as: 



2 x 2tt /* 72 sin 96,9 



COS t>T 



(5) 



47T 47T 

We must relate c to the observed ratio R, accounting for the 
anisotropy of disc and torus emission. Since the emission of ge- 
ometrically thin discs follows a cos 9 pattern, for a given viewing 
angle 9 V (calculated from the disc axis) the ratio between the real 
disc luminosity Lboi and the isotropic estimate L^oi is: 



L bo i 2 x 2ty JT' a cos # sin ( 

Liso 



4-7T COS 9 V 



2 cos 9 V 



(6) 



The ratio is smaller than unity for 9 V < 60°, thus for Type 1 AGN 
we likely have Lboi < Lboi- We are not able to determine cos 9 V 
for each source, but we can safely assert that ^ 9 V ^ 9t, since 
we are dealing with Type 1 AGN. Therefore a reasonable estimate 



cos# v ~ {cos6%-9 T 



Jo 



cos a sin f 



1 + cos 9t 



(7) 



A relation similar to Eq.|6]for the torus luminosity (Lt = Lir) is 
currently unknown. However, we can reasonably state that 



Lboi Lt ^ 

r iso r iso 
^bol 



(8) 



The lower limit corresponds to a thin disc-shaped emitting torus, 
the upper limit to an isotropic emitting torus. Both limits are 
rather unrealistic: the torus is expected to show a lower degree of 
anisotropy than the disc since we are able to detect radiation emit- 
ted from the side (i.e. Type 2 AGN); also, the torus is hardly an 
isotropic emitter since IR signatures are different in Type 1 and 
2 AGN (Calderone et al., in prep.). The above limits should then 
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Table 3. Results of our analysis. Columns are: (1) sample; (2) number of sources in the sample; (3) mean bolometric luminosity in the sample; (4) width of 
luminosity bin ( a : value in the first row is the standard deviation); (5) mean and standard deviation of LJ^ in the sample; (6) mean and standard deviation of 
parameter R = L^/Li® ,; (7) range of covering factor (Eq. |10| ; (8) range of torus opening angles; (9) range of Type 2 to Type 1 AGN count ratio (#2/#l). 
All means and standard deviations are computed using logarithmic values. 



Sample 


N src. 






logLjg 


R 


Cov. factor 


Q'j' 


#2/#l 


Whole 


3965 


45.72 


0.33 a 


45.18±0.27 


29+ ' 18 


0.54-0.70 


57-46 


1.2-2.3 


A 
B 
C 


408 
569 
389 


45.55 
45.80 
46.05 


0.10 
0.10 
0.14 


45.05±0.16 
45.22±0.16 
45.40±0.16 


31+ ' 14 
26+ 012 
22+ 10 


0.56-0.74 
0.51-0.66 
0.47-0.60 


56^12 
59^19 
62-53 


1.3-2.8 
1.0-1.9 
0.9-1.5 



bracket the real case. The amount of disc radiation intercepted (and 
re-processed) by the torus is: 



X 



7T-/2 



Lb 



Jo 



7T/2 



(9) 



cos a sin f 



Rearranging the previous equations, we find a relation between the 
observable parameter R = LJr /L™°\ and the covering factor c: 



1 + c 



<R<c 



(10) 



This relation can be inverted to find the allowed range of c and 6t, 
given a value of the observable parameter R. Finally, the covering 
factor c can be used to estimate the count ratio between Type 1 and 
Type 2 AGN: 

= Or = #2 #2 
#1 



(11) 



4ir #l + #2 #1 1-c 

The last three columns of Tab.[3]report the value of c, 9t and #2/#l 
corresponding to the observed values of R in all discussed samples. 



5 DISCUSSION AND CONCLUSIONS 

The main result of our work is the determination of the average 
covering factor of the torus using a very large data set. The ob- 
served fraction of IR to bolometric, isotropically equivalent optical 
luminosity is about 30%. This implies that the obscuring torus cov- 
ering factor is in the range 0.5-0.7 and that the opening angle f?x is 
40°- 60° . On average, our sources emit in the IR a similar fraction 
of their bolometric luminosities (~ 1/3). For each Type 1 AGN, 
there should be between 1 and 3 Type 2 sources. If there is a broad 
distribution of covering factors (as suggested by Elitzur, 2012) our 
Type 1 sample may be drawn preferentially from the lower end of 
the distribution. In this case our estimate of #2/#l ratio is a lower 
limit. The very basic prediction of the unified model that the torus 
re-processes a given amount of disc luminosity is verified (Fig. [TJ 
and Fig. O. The dispersion of this fraction is remarkably small, 
being at most a factor of 2. The broad-band spectral energy distri- 
bution (SED) from IR to near-UV are expected to be quite simi- 
lar among Type 1 AGN. A hint of the "receding torus" hypothesis 
is found in Fig.[TJ with logl^ oc 0.83 log L^. In Fig. [3] we 
show both data and model for a prototypical broad-band SED, in 
the three luminosity classes considered above (coded with colors). 
For each subsample we also compute a composite spectrum using 
SDSS spectra. 

At IR wavelengths the torus emission dominates. Spectral in- 
dices between the four WISE bands are very similar for different 
overall luminosities (Tab. |2). Despite the rather poor coverage, it 



appears that the IR emission is structured with at least two broad 
bumps. Such features are easily modeled by the superposition of 
two black bodies with temperatures of ~300 K and ~1500 K re- 
spectively. A naive interpretation is to consider the hotter one as 
emitted from the hot part of the torus facing the disc, at the dust sub- 
limation temperature. The colder one would come from the cooler 
outer side of the torus. This should be the region visible also in 
Type 2 AGN. 

The underlying optical continua are well described by a stan- 
dard Shakura & Sunyaev (1973) accretion disc spectrum. The 
dashed lines in Fig. [3] are the models of three accretion discs hav- 
ing the same bolometric luminosity as the spectra in the subsample, 
and masses 1.7 x 10 8 , 2.3 x 10 s and 3.4 x 10 8 M respectively, 
grossly in agreement with the (virial) masses calculated in S 1 1 . The 
WISE data points (in vL v ) lie a factor ~3 below the disc peaks (at 
log(^/Hz)~15.5). This factor corresponds to the value ~l/3 found 
in Tab. [3] and Fig. [2] The composite spectra follow closely the ac- 
cretion disc continuum in all but the lowest luminosity subsample, 
in which some other component is present at frequencies below 
log(i^/Hz)<14.9. This further component may be the starlight con- 
tribution from host galaxy (Vanden Berk et al. 201 1), as shown by 
the yellow line which is the sum of the accretion disc spectrum and 
an appropriately scaled template for an elliptical (quiescent) galaxy 
from Mannucci et al. (2001). At higher luminosity subsamples, the 
contribution from galaxy becomes relatively less important. 
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Figure 3. The disc-torus connection. AGN in our sample were divided in three subsamples according to their bolometric luminosity (see Tab.[3J, and associated 
to black, red and blue color respectively. The logarithmic mean and standard deviation of IR luminosities in each subsample are computed using data from 
WISE, and displayed as filled circles and error bars (points are slightly displaced for a clearer view). Composite optical/NUV spectra (solid color-coded lines) 
are computed as follows: SDSS spectra are transformed to rest-frame and de-reddened using Schlegel et al. (1998) and Pei (1992). The spectra are then 
rebinned to a common wavelength grid and a composite spectrum is computed as the geometric mean. The gray shades indicate the 68% level dispersion. 
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